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CONTRACTOR REPORT 


THE EFFECT OF A SMALL INITIAL CURVATURE ON THE FREE VIBRATION 

OF CLAMPED, RECTANGULAR PLATES 

SUMMARY 


An analytical method of obtaining the natural frequencies and mode shapes of 
clamped, rectangular plates having a small initial curvature is presented. Specifically, 
the perturbation technique is used to reduce the fourth-order plate vibration problem 
to the simpler membrane problem with modified boundary conditions that account for 
the bending effects. Excellent results are obtained when the bending rigidities are 
small compared with the in-plane forces. 


I. INTRODUCTION 


The problem of a plate having a small initial curvature has long been recognized 
[1], From small deflection theory, the total deflection is the sum of the initial deflec- 
tion and the deflection resulting from the applied load on a flat plate. However, few 
particular solutions are available in the literature. Timoshenko and Woinowski-Krieger 
[1] presented an exact analysis for the static deflection of simply- supported plates 
having a small initial curvature. The solutions to the vibration problems of such 
plates are, however, not available in the open literature. 

On the other hand, much work has been done on cylindrically-curved panels. 

The work of Sewall [2] revealed that natural frequencies of vibration of cylindrically- 
curved panels are obtainable in closed form for only two boundary conditions. Inter- 
estingly, it was found that the panel natural frequency deviates from the flat plate 
frequency only by a curvature term. This curvature term vanishes for large radii of 
curvature. However, Sewall’s formulation, which is complex and involves integrals, 
does not lend itself to easy application. Blevins [3] presented a simplified analysis 
based entirely on Sewall’s work. The Donnel equations for a cylindrical shell were 
used and mode shapes of a single- span beam satisfying pertinent boundary conditions 
were obtained. The Rayleigh-Ritz method was used to generalize Sewall's results to 
other boundary conditions. 

The problem presented here is different from those of Sewall and Blevins. We 
consider the free vibrations of a clamped, normally -prestressed plate with a small 
initial curvature using the method of singular perturbations. The governing equations 
are those for curved plates. There are no known analytical expressions for the free 
vibration of clamped, initially -flat plates even in the absence of normal prestress 
forces. Explicit expressions for eigenvalues and mode shapes for the free vibration 
of initially -curved, normally-prestressed plates have, however, been obtained in the 
present study, in the limit when the normalized bending rigidity is small. 

In the absence of shearing prestress, the asymptotic eigenvalue results for 
initially - curved , normally-prestressed plates are found to correspond to those of 
initially -flat, normally-prestressed plates. However, the eigenmodes are highly 
modified. It is shown here that the eigenmodes can be decomposed into the following 



three components: (i) the initial deflection, (ii) the flat plate deflection and (iii) the 

deflection resulting from the static equilibrium condition. 


II. PROBLEM FORMULATION 


We consider a thin rectangular plate having an initial deflection w o '(x* ,y'). If 

the plate is further deflected by an amount w^x'.y'.T), then the total deflection, 
w'(x' ,y',x), is 


w'(x',y',T) = w o '(x',y') + w^Cx’ ,y’ , t) 


( 2 . 1 ) 


The equilibrium equation given by Timoshenko and Woinowski-Krieger [1] is 
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for the case in which there is no shearing prestress. In non-dimensional form, equa- 
tion (2.2) becomes 
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where 


e 2 = D/N L 2 << 1 


(2.4) 


N q is the characteristic normal prestress value while L is the characteristic length 

used in normalizing the deflection w* and the co-ordinates x* and y*. Time is nor- 
malized using the characteristic frequency, u>, such that 


t = ojt and N o /yo) 2 L 2 = 1 


(2.5) 


where y is the mass per unit area of the thin plate. For a fully-clamped, rectangular 
plate, we require that w = 9w/9n = Oon all edges. 

If we let the dimensionless initial deflection w (x,y) be given by a double 
Fourier series of the form 
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w (x,y) = Ei E P m „ Sin(mrx) Sin (rmry/b) (2.6) 

o mxi 

m=l n=l 


where b is the dimensionless length of the rectangular plate in the y direction, the 
non-dimensional governing differential equation (2.3) becomes 
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where 


%n = -^ 2p mn /b2 > (b V + ^ 


( 2 . 8 ) 


It should be noted that the plate length in the x direction has been taken as the 
characteristic length , L . 

If we now make the substitution 


w t (x,y,t) = u(x,y,t) + <K x,y ) (2.9) 

in equation (2.7), we obtain the following dynamic equation: 
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with boundary conditions 


u(0,y,t) = u(l,y,t) = u(x,0,t) = u(x,b,t) = 0 
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u(x,y,0) = f(x,y) and (x,y,0) = g(x,y) 


(2.10) 
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and the static equilibrium equation: 
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with boundary conditions 


<K0,y) = >Kl,y) = Kx,0) = ^(x,b) = 0 


3 ifi 
9x 


(0 ( y) = f£(l.y) = fy - ( x >0) = f* (x,b) = 0 ) 


(2.13) 


Equation (2.10) is the governing differential equation for the free vibration of an 
initially- flat , anisotropically- prestressed rectangular plate, while equation (2.12) cor- 
responds to the static problem of a similar plate subjected to some generalized loading. 
Thus, the initial curvature effect is equivalent to that of some fictitious loading. 


III. METHOD OF SOLUTION 


The governing differential equations (2.10) and (2.12), subject to the condition 
equations (2.11) and (2.13), respectively, are solved by the method of singular per- 
turbations, in the limit when e -> 0 + . The procedure is to seek an approximate solu- 
tion valid in the core region of the domain. Thus, wa write 
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The resulting membrane-type governing differential equations, valid in the core 
region , are 
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Here, the subscript v denotes the order in e, while the superscript o denotes outer 
solution. Terms in which the dependent variable has a negative subscript are dis- 
carded. We note that the equations (3.2) for u^ 0 and are of second order as 
opposed to equations (2.1G) and (2.12) which are fourth order. 

It is easy to see that the expansions given in equation (3.1) cannot satisfy all 

the boundary conditions. Thus, the resulting solutions for u^ 0 and ifi °, which are 

equivalent to states of membrane deformation, will not be valid near the plate edges 
where bending deformations exist. Across these bending regions (boundary layers), 
the deformations change rapidly to satisfy the boundary conditions while the fourth 
order governing differential equations (2.10) and (2.12) are preserved. The solu- 
tions of these equations in the bending regions will then serve as modified boundary 
conditions for the membrane-type equations (3.2). 

For example, near x = 0, a stretching transformation takes the form 


X = x/e 


(3.3) 


In this case, we seek expansions of the 
u = u Q * + eUj* + e 2 Ug* + O (e^) 

4> = i|) Q * + e 4 >i + e 2 4>2 + 0 (e 3 ) 


form 
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The governing differential equations become 
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Equations similar to (3.5) can be written for the regions near x = 1, y = 0, and y = 
b. The dynamic and static problems will now be considered for arbitrary initial 
curvature . 
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3.1 The Dynamic Problem 


The dynamic problem for a fully- clamped , rectangular plate undergoing sinusoidal 
vibration can be described by the equation 


2„4 2 8 2 u 2 8 2 u 

e V u — & 2 7~2 


- A u = 0 
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(3.6) 


subject to 


u(0,y) = (0,y) = 0 


u(i,y) = d,y) = o ; 

0 < b < 1 

u(x,0) = (x,0) = 0 ; 

u(x,b) = |-£ (x,b) = 0 . 

3y J 


(3.7) 


It is easily verified from the literature [4,5,6] or otherwise that an approximate solu- 
tion to equation (3.6), subject to the boundary conditions (3.7), can be written as 


u(x,y) = A q sin (mrx) sin (miry/b) + e A q {(mr/ 3 1 )(2x- 1) cos (nirx) sin (rmry/b) 

+ (mir/ 3 2 b 2 )(2y-b) sin (mrx) cos (rmry/b)} + e 2 2 A q {(n 2 n 2 / 6^ 2 )x( 1-x) 

cos (mrx) cos (rmry/b) + (n-n723 1 2 3 2 )(23 2 b- 6p( 2x “l) cos (mrx) sin (mrry/b) 
+ (mrr/26 1 3 2 2 b 3 )(23 1 - 6 2 b)(2y-b) sin (mrx) cos (miTy/b) 

+ (m 2 ir 2 /62 2 b 4 ) y(b-y) sin (nirx) sin (miry/b)} + e {(z^/g^) exp (-gjX/e) 

+ (z 2 / 3 X ) exp (-3 1 (l-x)/e) + (z 3 /3 2 ) exp (-3 2 y/e) 

+ (z 4 /3 2 ) exp (”3 2 (b-y)/e)} + ^ ex P (" 3jX/e) 

+ (Zg/Pj^) exp (-3 1 (l-x)/e) + (z ? /3 2 ) exp (~3 2 y /e) 

(3.8) 

+ (Zg/3 2 ) ex P (-6 2 (b-y)/e)} + O (e 3 ) 
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while the eigenvalues X 


2 


are given as 


X 2 = Bj 2 n 2 2 + B 2 2 m 2 ^ 2 /b 2 + e 4 {$ 2 n 2 tt 2 + B 2 2 m 2 Tr 2 /b 2 } 

+ e 2 {(7r 4 /b 4 )(n 2 b 2 + m 2 ) 2 + (12 7r 2 /b 4 ) (n 2 b 4 + m 2 )} + O (e 3 ) . (3.10) 


The eigenvalues for plates with inverse aspect ratios, b, varying between 0.1 
and 1.0 and for the dimensionless normal prestress values Bp B 2 between 0.1 and 

1.0 are presented in Figures 3.1 through 3.4 for values of e ranging between 0.0010 
and 0.2500. 

It is seen from Figure 3.1(a) that the zeroth and second order solutions for 
b = 1 (square plate) are essentially indistinguishable for values of e of 0.0025 and 
0.0010. However, as b decreases, a considerable difference develops between these 
two curves, indicating the solution to be more sensitive to the value of the normalized 
bending rigidity, e. A similar trend, albeit more pronounced, is observed for 
decreasing values of the normalized prestress Bj^ and B 2 as shown in Figure 3.1(a)iii 

and Figure 3.1(b). It is observable from Figure 3.1(b)iii that the difference between 
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2 

the second order and zeroth order eigenfrequency ( X ) can be as much as four orders 
of magnitude at mode numbers m = n = 15 for values of e = 0.0250 and higher. 

Figure 3.2 shows the second order eigenfrequencies for e = 0.0010 plotted 
against the mode numbers with the inverse aspect ratio b as parameter. The observed 
trend of decreasing values of the eigenfrequency with increasing b is easily inferred 
from the eigenfrequency expression equation (3.10). 

The same trend of decreasing eigenfrequency with increasing b is also evident 
in Figures 3.3(a) and (b) in which the additional effects of combinations of m and n, 
the mode numbers in the x and y directions respectively. Figure 3.3(a), and com- 
binations of and gg* the normalized prestress in the x and y directions respec- 
tively, Figure 3.3(b), have been highlighted. 

Finally, Figures 3.4(a) and (b) show the eigenfrequencies plotted against the 
normalized prestress for various combinations of values of m, n, and b. At higher 
values of m and n, Figure 3.4(b), the eigenfrequency is essentially independent of 
6^ and $2 at the higher e values. 

Next, we consider the static problem for an arbitrary initial curvature. 


3.2 The Static Problem 


The governing differential equation for a fully-clamped, rectangular plate takes 
the form 


2 U 4. o 2 3> 
e V ip - g 1 — \ - 
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subject to 


<K0,y) = |f (0,y) = 0 
<Ki,y) = d,y) = o 
H x , 0) = (x , o) = o 

^(x ,b) = |f (x,b) = 0 


0 < b < 1 > 




(3.11) 


(3.12) 


Using the particular form of w Q (x,y) given in equation (2.6), and by compari- 
son with equation (2.12), the right hand side of equation (3.11) is seen to take the 
form 
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q mn Sin (nirx) Sin (muy/b) 


where q mn was defined in equation (2.8) in terms of P mn > the Fourier coefficients 

obtained by expressing any function f(x,y), describing the initial curvature, as a 
double Fourier series. The procedure is to solve equation (3.2b) for v = 0, 1, 2,... 
subject to appropriate boundary conditions. It is easy to see that the boundary value 
problem for \p Q takes the form 


e 


2 

1 





q mn Sin (nitx) Sin (miry/b) 


(3.13) 


subject to 


<l^ o (0,y) = i|> 0 (l,y) = <|< o (x,0) = i|) Q (x,b) = 0 . (3.14) 


Clearly, the solution to equation (3.13) satisfying all the conditions of equation 
(3.14) is 


4> 0 = Y mn sin (nmc) Sin (miry/b) 


(3.15) 


with 


,, 2 2 2 „ 2 2 2 ,. 2 . 

Y mn = <W (e l " it + g 2 m ir /b ) 


(3.16) 


To the next order of approximation, ^ must satisfy the differential equation 
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tjd.y) = (-l) n ^ 1 (0,y) 

^(x,0) = - y mn (mir/b) Sin (nirx) 


(3.18) 


^(x,b) = (-l) m ^(x.O) 
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In previous studies [7,8], equation (3.17) subject to equations (3.18) was solved 
numerically. In general, any of the standard numerical methods of solving parabolic 
partial differential equations such as the Gauss- Seidel, Jacobi and the under- and 
over- relaxation (SOR) methods can be used. However, using the singular perturba- 
tion analysis technique, analytical expressions can be obtained to leading and higher 
orders with increasing accuracy as e + 0 + . It can be shown that 


\p^ = {A^ cosh (a^y) + A 2 sinh (a,y)} Sin (mrx) 

+ (A 3 cosh ( c* 2 x ) + A 4 sinh ( a 2 x )) s * n ( m7T y/b) 


where 


(3.19) 


a 1 = n tt 6 ., / e 2 and a 2 = mT7 8 2 ^l b 


(3.20) 


and the constants A^, ..., A^ are observed to take the form 

A 1 = ' mTr W 3 2 b 

A 2 = - A^jt- l) m + cosh ( ct.|b) }/sinh (ct-jb) 


> 


A 3 = ” n7T Y mn /3 1 

a 4 = - A 3 {(-l) n + cosh (a 2 )}/sinh (a 2 ) 


(3.21) 


The method of solution is further extended to The boundary value prob- 

lem to be solved for is 


e 


2 

1 



+ 



(3.22) 


subject to 

P 2 (0,y) = - (A 1 cosh (o^y) + A 2 sinh (o^yjjniT/^ - (A 4 c^/^) sin (rmry/b) 
^ 2 (l,y) = (A 1 cosh (a 4 y) + A 2 sinh (c^y)) (-l) n nir/gj 

+ (a 2 /8 4 ) (A 3 sinh (a 2 ) + A 4 cosh (a 2 )) sin (m-ny/b) 
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i|* 2 (x,0) = - (A 2 a 2 /£ 2 ) sin (n-rrx) - (rmT/bB 2 ) (A 3 cosh (a 2 x) + A 4 sinh (a 2 x)) 

^ 2 (x,b) = (A 3 cosh (a 2 x) + A^ sin (a 2 x)) (-l) m mn/g^b 

+ (a 4 /3 2 ) (A^ sinh (otjb) + A 2 cosh (a^b)) sin (mrx) . (3.23) 

We seek the solution to equation (3.22), subject to equations (3.23), in the 

form 


*2 = ha + hb 

where 

^2 a = cosh ( a i3p + &2 sinh (a l y)] sin ^ n7TX ) 

+ [g 3 cosh (a 2 x) + g^ sinh (a 2 x)] sin (rmry/b) 

+ y sin (n-rrx) sin (mny/b) 

with 

Y = - (nV + m 2 rr 2 /b 2 )/(B 1 2 nV + 8 2 2 mV/b 2 ) 

and 

^1 A 2 a 1^2 

g 2 = [A^ sinh (ctjb) + 2A 2 cosh (a-jb)]/^ sinh (a^b) 

g 3 = " A 4 a 2^1 

g 4 = a 2 [A 3 sinh (ag) + 2A 4 eosh (a 2 )]/P 1 sinh (a 2 ) 
Also, the expression for takes the form 


(3.24) 


(3.25) 


(3.26) 


11 


^ 2 b = [E 1 Sinh (a l y) + E 2 COSh (a l y) + E 3 y sinh (a l y) + E 4 y cosh (a i y ^ sin ( n7rx ) 

+ [Eg sinh (a 2 x) + E g cosh (eye) + E ? x sinh (a 2 x) + E g x cosh (ye)] sin (rrnry/b) 

(3.27) 

where 


E^ = - E 2 [ ( - l) m + cosh (a 1 b)]/sinh (a^b) - E g b - E^b coth (a^b) 
E 2 = - mir [Ag K 1 + A 4 Kg] / 3 2 b 
E g = A 4 3 4 n 2 tt 2 / 3 2 2 a i 

E 4 = E 3 A 2 /A l 

Eg = - Eg [ ( - l) n + cosh (a 2 )]/sinh (a 2 ) - E ? - E g coth (a 2 ) 

Eg = - rnr [Aj K g + Ag K 4 ] / 3 1 
E ? = Ag 3 2 m 2 tt 2 /3 1 2 a 2 

E 8 = E 7 A 4 /A 3 






J 


while 

n + 1 9 9 9 

= nTr[ 1 + (-1) cosh (a 2 )]/(a 2 + n it ) 

Kg = (-l) n+1 nu sinh (a 2 )/(a 2 2 + n 2 tt 2 ) 

m+1 9 9 9 9 

Kg — m tt [ 1 + (-l) n cosh (a 1 b)]/b(a 1 ^ + it / b 9 

K^ = (-l) m+1 [Bit sinh (a-jbVMa^ 2 + m 2 Tr 2 /b 2 ) 

9 

This completes the solution of the problem to O ( e ) . 
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(3.28) 


(3.29) 
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IV. CONCLUSION 


The small initial curvature effect on the free vibration of clamped, rectangular 
plates has been presented., The eigenfrequeneies for plates with inverse aspect ratios 
varying between 0.1 and l.'O and for the dimensionless normal prestress between 0.1 
and 1.0 have been presented for values of e ranging between 0.0010 and 0.2500 in 
Figures 3.1 through 3.4. 

It is established here that a small initial curvature has no effect on the fre- 
quency of vibration of the plate. However, its effect is manifested in the eigenmodes. 
The eigenmodes obtained are of increasing accuracy as e 0 + , similar to the trend 
for the eigenfrequeneies presented in Figures 3.1 through 3.4. 

The solutions obtained from above provide a guide for researchers that imple- 
ment the finite difference, finite element or other numerical schemes on the computer 
for the solution of curved plate problems. This utility of the singular perturbation 
technique in providing simplified expressions which yield solutions valid in the 
asymptotic limit and which can serve as useful guides for the pure numerical analysis 
approach has been attested to by Nayfeh and Kamat [7] and Ramkumar et al. [8]. 
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